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. 1. Introduction 

The seed of the phase-space formulation of Quantum Mechanics dates back to 1932, 
when Wigner introduced his famous quasi-probabihty function in phase space [Ij: 

The relation between Wigner's function and Weyl's correspondence was understood 
by Groenewold [2] and fully developed by Moyal |3j, who established an independent 
formulation of Quantum Mechanics in phase space for the canonical Poisson bracket 
{q,p} = 1. In this equivalent, 'statistical' or hydrodynamical formulation of Quantum 
Mechanics, the expectation value of an operator is computed as the average of the 
corresponding function on the phase space with the 'probability' density given by the 
Wigner function. 

In light of these results, it was suggestive to look for an autonomous formulation 
of Quantum Mechanics in phase space involving only classical functions and valid for 
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any Poisson bracket. This is the so-called deformation quantization program proposed 
by M. Flato and collaborators in the 1970s The problem is stated as follows: given 
a smooth manifold M with Poisson bracket {-,■}: C°°(M) x C°°(M) C°°{M), find 
an ]R[[/i]] -bilinear associative deformation of the point-wise product on C°°{M), -k, such 
that 

f^9-9^f = Hf,9} + Oiff), V/, g G C^iM). 

The existence of i^r-products for any Poisson manifold was a longstanding problem 
in Mathematics solved in 1997 by M. Kontsevich [5J. As a consequence of his more 
general Formality Theorem, he showed that any Poisson manifold can be quantized by 
deformation and classified the ^-products (see globalization aspects in |6]). 

One of the main advantages of the deformation quantization prescription (or 
equivalently, the phase-space formulation of Quantum Mechanics) resides in the 
quantum-classical transition problem. The classical limit is obtained in a clear and 
mathematically rigorous way by taking h —>■ 0. 

The phase-space formalism has been applied successfully to the description of a 
spinless particle ([71 E]), whose phase space is with the canonical Poisson bracket. 
More complicated is the phase-space description of a spin, the phase space being the 
sphere S'^ . The phase space of a spin was given in terms of the atomic (or spin) coherent 
states, O [ini [m [I2] , which in a sense generalizes the Cahill-Glauber construction for 
the spinless case [T3] . 

Although the results of Kontsevich guarantee the existence of a T*r-product for 
any classical system (i.e., for any Poisson manifold), his explicit formula is rather 
complicated and, for example, it is not obvious how to take advantage of the symmetries 
of the system. Varilly and Gracia-Bondia (see [H] for SU{2) and [15] for any compact 
group) showed that the appropriate setup for systems with symmetries is given by the 
Stratonovich-Weyl (SW) correspondence [H] (see Section [2]). In addition, they proved 
that this approach is equivalent to the spin coherent state representation. These works 
were generalized by Brif and Mann JlTl [18], establishing the SW correspondence for 
systems with an arbitrary finite-dimensional Lie-group symmetry. 

In the present work we tackle the problem of the dynamics of spin systems in 
the phase-space formulation and its classical limit. Recently Klimov and Espinoza [T9] 
derived a differential form for the star product in the spin case. On the other hand, 
in the spinless case the evaluation of the star product becomes simpler with the use 
of Bopp operators [20]. We use Klimov and Espinoza's results and work out the Bopp 
operators for spins. Takahashi and Shibata already found this generalization for the 
specific cases of normal and anti- normal ordering [TOj [11]. Herein, we re-obtain their 
results and generalize them for an arbitrary ordering. In particular, we obtain the Bopp 
operators for the important case of symmetric ordering. We will apply our results to 
show how the use of Bopp operators simplifies the derivation of the dynamical equations 
for spins. 

Finally, we deal with the problem of a system in contact with a bosonizable bath 
as a model of quantum dissipation [2]]|. Classically, the effective dynamics of the 
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system is described by Langevin or Fokker-Planck equations [22] • In the quantum 
domain and under certain conditions (essentially weak coupling between system and 
environment), the dynamics of the system can be formulated in terms of a quantum 
master equation for the (reduced) density matrix [23]. For spinless particles the phase 
space transform of the master equations yields the quantum generalization of the Klein- 
Kramers equations [21]. In the sense of Caldeira and Leggett, who pose the open system 
dynamics as a quantization problem [25] , the phase-space formulation "closes the circle" , 
giving a quantum version of the Fokker-Planck equations. With the help of the Bopp 
operators for spins we easily obtain quantum Fokker-Planck equations, recovering the 
corresponding classical Fokker-Planck equations for the rotational brownian motion [26] 
in the limit h —>■ 0. This provides a natural framework to link the classical and quantum 
theories of dissipation. 

The paper is organized as follows: 

Section [2] is a brief survey on the phase-space formulation of Quantum Mechanics 
from the point of view of the Stratonovich-Weyl postulates. The general results are 
illustrated by the example of the non-relativistic spinless particle, introducing the 
definition of Bopp operators. 

Section [3] deals with the phase-space formalism of spin systems. We generalize 
the Bopp operators and tackle the problem of the quantum-classical transition in 
these systems. In addition, we work out in detail the particular case of quadratic 
Hamiltonians. 

In Section H] we apply our results to open quantum spin systems. We transform 
into phase space the density matrix equation for a spin in contact with a thermal bath, 
obtaining quantum Fokker-Planck equations. This transformation becomes simple with 
the help of Bopp operators. We write explicitly the dissipative equations for the linear 
Hamiltonian [isotropic spin) and compare it with its spinless analogue, the damped 
harmonic oscillator. Finally, we take the classical limit for a general quantum Fokker- 
Planck equation, recovering the corresponding classical equation. 

2. Phase-space Quantum Mechanics 

In this section we review some basic facts and results on phase-space quantization. The 
material is now standard and we closely follow the conventions and notation of [T7t [I8]. 

2.1. The SW postulates 

Consider a physical system possessing a group of symmetries described by a finite- 
dimensional, connected and simply connected Lie group G. We denote the phase 
space of the system by V and assume that G acts transitively on V. That is, for 
any Xi, X2 G V, there exists g E G such that g ■ xi = X2. Let H be the Hilbert space of 
our system, 0(TC) the set of operators on H and U : G ^ 0(H) an irreducible unitary 
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representation of G. The Stratonovich-Weyl correspondence ([I5|) is a cr-parameterized 
map W^"'^ : OiH) C°°{V) satisfying the following properties for any A, B E 0{1-L): 

(i) Linearity: 

A f— > W^'' is linear and bijective. 

(ii) Reality: 

^At^(^) = Wa\^Y-> where A'^ is the adjoint of A. 

(iii) Standardization: 

I W^^\x)d^i{x) = Tr^. 
Jv 

(iv) Traciality: 

/ W^^\x)w''s''\x)d^i{x) = Ti{AB). 
Jv 

(v) Covariance: 

W^/J{x) = W^/\g-x), ygeG. 

where /i is an invariant measure and A^ := hl{g^^)AU{g). 

As proved by Wigner et al. [7j for spinless particles, by Varilly and Gracia-Bondi'a 
for spins [H] and recently by Brif and Mann pT| IT8] for systems with an arbitrary 
Lie-group symmetry, properties (i)-(v) uniquely determine the SW correspondence. 

Reality condition (ii) ensures that W^'' is real whenever A is hermitian and property 
(iii) is simply a normalization. Linearity and Traciality conditions are essential to the 
formal interpretation of Quantum Mechanics as a statistical theory. Let w'f'* be the 
image by the SW correspondence of the density matrix of the system, g. For any 
operator A, {A) = Tt^Aq) and due to (iv) we obtain: 

(A) = [ W^/\x)W^g~''\x)dnix). (1) 
Jv 

Therefore, the expectation value of an observable A can be computed as the average on 
V of W^^ weighted by Wg~"\ However, Wg'^^ may be negative at some points ([7j), so 
that in general it does not define a true probability distribution. For this reason 
is sometimes called a quasiprobability distribution. 

Finally, the covariance property (v) means that the SW correspondence commutes 
with the action of G, thus preserving the symmetry. 

Notice that under the SW correspondence, an operator is mapped to a set of 
functions labeled by cr, which is related to the operator ordering prescription. In 
particular, for a = 0, 1 and —1 are known, respectively, as symmetric, normal 

and antinormal functions. The nomenclature is due to the fact that, for the spinless 
particle, cr = corresponds to the symmetric ordering (|3J), whereas a = 1,-1 are 
associated to the normal and antinormal orderings, respectively ([27]). 
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2.2. The SW kernel 

The vector space 0(TC) equipped with the bihnear form {A, B) = Tr(A^i?) is a Hilbert 
space and due to the Riesz theorem, there exists a cr-parameterized set of operator- 
valued functions on V, A^'") G C°°{P) (g) 0{n), called the SW kernels, such that: 

W^/\x) = Tt{AA^''\x)). (2) 

The properties (ii)-(v) introduced above are equivalent to the following conditions 
on the kernels ([12 [13 [H]): 



fii') 



fiii') 



(iV) 



fvO 



A(")(2;)d/i(x) = 1 
A^^^g ■ x) = A^''\xy~\ WgeG. 



Remarkably, the tracial property (iv') allows to invert the SW map ([2]), yielding a 
generalization of the Weyl rule: 

A = [ w[^\x)A^-''\x)dfx{x). (3) 
Jv 

Notice that the kernel involved in ([2]) is A^'^^ whereas the kernel entering formula ([3]) 
IS A^"'^). The case CT = is privileged and A'-'^^ is said to be a self-dual kernel. 

2. 3. The star product 

Once the Stratonovich-Weyl correspondence has been constructed, the connection with 
the notion of i^r-product is straightforward. The T»r-product is obtained by transferring to 
V the associative algebra structure of C(7i) through the SW map. That is, 

{W^/Kwi^^){x):=W^/^{x). (4) 

for any two operators A and i^. 

At this point, a natural question arises. Given an arbitrary ^-product on C°°{V), 
is there a Poisson bracket associated to it in a canonical way? The answer is positive 
and the proof is easy. Using associativity, a direct computation shows that 

{f,9} = ^^^^^ mod /I, f,geC^iV) (5) 

I A more general relation mixing different orderings can be given (see for example [T5J), but we will 
not discuss it in this work. 
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is a Poisson bracket. This is a nice way to see that any quantum system has a 
corresponding classical system which is obtained in the limit h ^ hj replacing 
commutators by Poisson brackets. For example, we can write the transformation of 
the von Neumann equation to phase space: 

where H is the Hamiltonian of the system. Hence, the star product determines the 
dynamics on the phase space and reproduces the classical Hamilton equations in the 
limit h 0. 



2.4- An example: the spinless non-relativistic particle 

The classical phase space of the spinless non-relativistic particle is with canonical 
Poisson bracket in coordinates 

{g,p} = i- (7) 

The Hilbert space of the quantum system is L^(M) and its dynamical symmetry 
group is the Heisenberg-Weyl group i/3. The Lie algebra of is generated by three 
elements /, a and with Lie brackets: 

[a, /] = [a^ /] = 0, [a,a^] = /. (8) 

i.e. and a are creation and annihilation bosonic operators, respectively. 

The SW correspondence in this case maps each operator p) into C°°(M^), where 

dib 

(#)(go) := go^(go), {p^){qo) ■= "i^^l^o, Vgo e M (9) 
for any if; G ^^(IR). 

Using the standard complex coordinate z = {q + ip)/ V2h the SW kernel reads: 

AM(z) = I/tt / etl«l'e«*^-«^*e«'^^-«*M2e. (10) 
Jc 

Let us focus on the case a = (symmetric case). For this particular value of the 
parameter a, ffTOl) becomes the famous Wigner function [1]: 

/+00 
eWft(g_iM|A|g+iM)dM. (11) 
-00 

The T*r-product is the Moyal product, whose closed form was introduced by 
Groenewold in the forties [2]: 

<i(.p)^<'e— r.||-J|. (12, 

A useful representation of the star product was introduced by Bopp [20]. He found 
that, for any two operators A{q,p) and B{q,p), the i^r-product can be casted in the 
following form: 

<^<) = A(Q,P)<) 

^ ih d ^ ih d , ^ 
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where now A{Q,P) should be understood as an operator acting on C°°{V). The 
operators Q and P are usually known as Bopp operators. Observe that [Q, -P] = i^, 
which can be checked directly from ( |T3l) . As a consequence, Q and P are operators on 
C°°{V) satisfying the canonical commutation relations. 

Using the reality condition of the SW kernel we can write the time evolution of the 
symbol of the density matrix, Wg'^\ in a simple way: 

= ^Im(iJ(Q, P))iyf (14) 

with H the Hamiltonian of the system. The classical limit is very easily obtained now. 
Expanding up to first order in h, 

HiQ, P) = H{p, q) + '^{d,H^^ - d,H^) + 0{h') (15) 

we recover the classical Liouville equation: 

= {if,iyf }. (16) 



3. The spin case 

We now turn to spin systems, our main interest in this work. The dynamical symmetry 
group is SU(2), and the phase space is the unit sphere S^. The classical Poisson bracket 
is suitably written in terms of the components of 

m(6', 0) = (sin 6^ cos 0, sin 6* sin 0, cos 6*) . (17) 

Nameljli, 

{murrij} = ^eijkrrik, (18) 

where eijk is the Levi-Civita symbol. Strictly, the Poisson bracket should be {mi,mj} = 

geijkrrik, where g is the gyromagnetic ratio, which is related to the quantum spin number 

by S = fis/gK with /i^ is the magnetic moment. From now on we set /i = /i^ = 1, so 

that S = 1/g. The classical Liouville equation may be expressed ([2B]) as 

If)/ \ f)J-f 

dtW^/^ = {H, W^;^} = --^ ■ (m X Beff Beff := , (19) 

^ ^ b am V / ^ am 

with the divergence {d/dm) ■ A = "^-dAi/dmi. The purpose of this section is to 

generalize the Bopp operators for the spin case, which simplifies the study of the 

dynamics in phase-space reported in the literature [28] . 

3.1. The SW kernel and the star product 

The SW kernel reads [El [H [HI [H] : 

I — I ' 

A('^)(^,0) = J^^5^(5,5;/,0|5,5)-- J2 TimY;je,<p) (20) 

1=0 m=—l 

§ Summation over repeated indices must be understood throughout the paper. 



Bopp operators and phase-space spin dynamics 



8 



where (j, j; 0|j, j) are Clebsch-Gordan coefficients, Ti^ = a/ (2/ + 1)/ (25* + 1) x 
^^■^.,(5', j; 5, j'lS*, j')!^, j')(5', j| are the irreducible tensor operators and YJ^ are the 
spherical harmonics. 

Recently, Klimov and Spinoza derived a differential form for the star product in 
the spin case [19]: 



where is the Casimir operator on the sphere, 

1 



de^^'''' de~' sin^ 6 86^ 



C^Yi^ = l{l + l)n 



F(e)Y,„, = F{l)Y,„„ 



k=0 



86 sin 6 



F{1) = ^/{2S + 1 + 1)1(23-1)1, 

a.= i-^y 

j\{2S + j + l)V 



and Ns = V2S + 1F^{0). 



(22) 
(23) 
(24) 



3.2. Bopp operators in the spin case 

The evaluation of ( l2TI) seems quite involved. Here, we show that it can be simplified by 
generalizing the Bopp operators ( fT3l) for the spin case. That is, we look for operators 
S^'^^ acting on C°°{V) such that: 

z = 1,2,3. (25) 



A ' 



Let us invoke the expression derived in the Appendix for -^Yi^, Eq. (]A.13p . 
Noticing that = cos6, L3 = —id^, i(m x L)3 = sindde and recalling that every 
function on the sphere is a linear combination of spherical harmonics we deduce that 
for any operator A: 



m3r^i^)(£^) + i(mxL)3r/r(/:^) + ^L3 



S3 " A 

where f,f\C^)Yim = V^"\l)Yi^, i = 1,2 and 
F^-''{1) 



W 



(^) 



(26) 



(0 



2(2/ + 1) 
F^--{1) 



F''~\l + 1) (25 + / + 2) (/ + 1) - F^-^l - 1) (/ - 25 - 1) / 
F"-^(Z + 1) (25 + / + 2) + F"-^(/ - 1) (/ - 25 - 1 



(27) 



2(2/ + 1) 

The vector m has been defined in ( ITTl) and L is the angular momentum operator acting 
on the sphere. Explicitly, 

d 



— 1 m X 



dm 



which satisfies [Li,Lj\ = ie^fc-^fc- The following identities are useful: 
[rrii, Lj] = ieijkirik; [(m x L)^, Lj] = ieijk{ni x L)^ . 



28 



(29) 
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5f ) = msvnC') + i(m x L),fjri^") + ^L,. 
The covariance of the star product at the infinitesimal level implies that 



1 



whence 



Using Eq. (ED) along with (EHD and ((291) we obtain 



1 



r^r(/:^) + i(mxL),77r(^') + 2^. 



(30) 
(31) 
(32) 

(33) 



The operators S^'^^ are the sought generalization of the Bopp operators to the spin 
case. Now, by the associativity of the star product, we can extend to any function 
^(5*1, S2, S3) of the spin operators: 

W^i^i = A(5f\5^^),5f))W^('^) (34) 

where, equivalently to the spinless particle case (see Eq. (HI), v4(5f\5f\4''^) is the 
same function of its arguments as ^(5*1, 5*2, S3), now acting on C°°{V). In order to ease 
the notation we will simply write A{S^'^^) instead of A{S["\ s!f\ 33°^^) . 

The von Neumann dynamics is governed by the equation (compare with (fT4l) ): 

= 2Imfif(5("^))iyf ) (35) 



(^) cWi 



Besides 



In this sense Bopp operators are the phase-space 



analogues of the corresponding Hilbert-space operators. 
Remark: 

In the a = 0, ±1 cases we can write: 

S 



(^) 
V2 



1 



2(2/+l) 



(/ + l)v/(25 + l)2-(/ + l)2 

^ + 1 



/v/(25 + l)2 + /2 



a = -I 

(7 = 

a = 1 



1 



2(2/+l) 



_ 1 

2 

v/(25 + 1)2 - (/ + 1)2 - ^(2^ + 1)2 + P 
1 

2 



a = -1 
a = 
a = 1 



(36a) 



(366) 



The form of S^"^^^ was already found by Takahashi and Shibata in [1 Ij by extending the 
construction of Cahill and Glauber. We have recovered their result and generalized it 
to any value of a in terms of the symbol ([2]) for the differential form of the star product 
of Klimov and Espinoza (12T|) . This covers the important case of symmetric ordering. 
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3.3. Large S 

We turn to the large S (or classical) limit. Firstly, notice that for 5 ^ 1 we can write 

m-. 



Inserting the above expansion in fl27l) we inmediately find that: 



(37) 



2S 



a + 



1 



(a^ - 1) + O 



1 



(38) 



2(25 + 1)' ' \S\ 

The above expressions allow us to evaluate the quasiclassical evolution of Wg"^^ in a 
simple way (see also [29]). This is an alternative to the spin coherent path integral 
formalism used in [30], [31]. In the classical limit, fj['^^ = S + {a + l)/2 + 0(1/5') and 

^2 



^^"^ — cr/2 + 0(1/5). Let H{Si, S2, S^) be the Hamiltonian. Also in the classical limit. 



lm{H{S 



-i dH 
25 9^ 



and using that for any b G M"^, 
-ib-L ^ 



dm 



m 



X b]. 



(39) 



(40) 



we recover the classical Poisson bracket (HM . obviously independent of a. 



3.4. Quadratic Hamiltonians 

We apply now the formalism to the following case: 
H = DjiSiSi BjSi, Di 



- D,, (41) 

which is the most general Hamiltonian quadratic in the variables 5j. This 
example covers many applications in magnetism [32] and quantum optics (see, 
e.g., [331 EH). With the help of ([31 and ^ we obtain that Im(if) = 



— 1 



DjkLj [mkr]['^^ + i(m x L)fc?72'^'' j + BjLj/2 and using fHU]) we finally get: 

1 d 



S dm 



m X B 



cflf 



s 



with 



Q 1 

eSij 



[Beff],- + 2SD,k mfc(77f ) - 5) + i(m x L)^//^' 



(42) 



(43) 



Equation 



is a compact way of writing the quantum dynamics in phase-space. 
B|^ is defined in such a way that the first term gives the classical effective field 
[Befr]j = SBj + 2S'^Djkmk (see Eq. f lT9l) ) and the second one the quantum contributions. 
When Djk = the quantum time evolution equation ( H2l) coincides with the classical 
one ([m Theorem 5]). This is analogous to the harmonic oscillator for the spinless 
particle [7] where the classical and quantum equation of motion in phase space are 
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the same. Notice, however, that ahhough the phase-space time evolution equations in 
the classical and quantum cases are formally identical, the allowed solutions for the 
differential equations are different at the classical and quantum levels. 



with r]['^ and given in (1271) . For cr = ± 1 we recover the results found by Klimov (see 
equations (65) and (67) of Ref. [2S]). However, we obtain a different result for cr = 
(Eq. (66) in Ref. [28]). 

Finally, using the asymptotic forms of fj[^^ and fj^2^\ Eq. (1381) . we obtain the 
quasiclassical form of fH2l) . In the case H = —DS^ the result agrees with derivations 
previously appeared in the literature (see the Conclusions of Ref. [29]). 

4. Application to Quantum Master Equations 

The field of open quantum systems deals with the unavoidable interaction between 
the system of interest and the environment pi]. Under assumptions such as weak 
coupling between the environment and the system and fast bath dynamics [35l [25] , the 
environment may be modelled as a collection of oscillators. With these hypotheses one 
can write the total Hamiltonian (system plus environment) as: 



where a is an oscillator index and the coupling terms F are functions of the system 
variables. 

Classically the dynamics of the system degrees of freedom is formulated in terms 
of Langevin or Fokker- Planck equations [22l[36]. At the quantum level, Qa and Pa are 
the position and momentum operators and F is a hermitian operator which depends on 
the system variables. The total system is unlikely to be in a pure state and a density- 
matrix description is required. For observables depending only on the system variables, 
the required object is the reduced density operator g = Trbath(ftot), where one traces the 
bath variables out. For weak system-bath coupling a closed dynamical equation for g 
can be obtained by perturbation theory. This is the case of many problems in quantum 
optics, chemical physics or magnetism [21]. As far as the classical limit is concerned, a 
semi classical quantum master equation would be sufficient. Let T be the temperature of 
the bath and 7 the damping coefficient (measuring the coupling strength). Considering 
that 7/(5T) < 1 one arrives to [21123]: 



Remark: 



Consider the particular case H = —DS^. Then, Eq. ( l42l) simply reads: 




(44) 




(45) 




(46) 
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where h.c. means "hermitian conjugate". The term —i[Hs,g\ gives the isolated-spin 
unitary evolution and the rest encodes the bath influence. 

Now, we address the problem of working out the phase-space transform of ( l46l) . In 
the spinless case for the generic Hamiltonian Hg = p^ /2m + V{q) and bilinear coupling, 
F = q, the master equation (H6l) is the celebrated Caldeira-Leggett equation [21]. Then 
the phase-space transformation gives the first quantum corrections to the Klein- Kramers 
equation, which is nothing but the Fokker-Planck equation for a particle in a potential 
[36]. 

In order to obtain the quantum corrections to the classical Fokker-Planck equations 
for spins, we will make use of the generalized Bopp operators. Using the results of Section 
13.21 for transforming Hilbert space operators into operators acting on C°°{V) we obtain: 

2 



and 



[F, Fg] + h.c. 



[F,[H,,F]g]+h.c. 



c.c. 



g 



(47) 



F{S 



where c.c. means "complex conjugate'' 
formulation: 



c.c.) ([i7s(5('^)),F(5('^^)] - c.c.)!^^'^) (48) 
Then, Eq. (HGI) reads in the phase-space 



I- g 



2Im(i/,(5('^))) + 



(49) 



+ 47T 



Im^FiS^) - —lmiFiSn)lm{[HgiS^'^^),FiSn]) 



This is the phase-space equivalent of the master equation ( 1461) and should be viewed 
as the quantum generalization of the classical Fokker-Planck equations. We emphasize 
that the transformation to the phase-space formalism has been rather simple thanks to 

(El. 



4.1. The open isotropic spin versus the open harmonic oscillator 

To gain some insight we particularize Eq. (1491) to an specific example. Our intention 
here is to highlight the differences between the spinless particle and the spin system. 
For that, we closely follow the work of Caldeira and Leggett on the quantum master 
equation for a spinless particle [24j. In this case the phase-space transform of the last 
two terms of (l46l) gives the dissipation and diffusion terms of the classical Fokker-Planck 
equation [25] . In particular, for the harmonic oscillator, both the Caldeira-Leggett and 
Fokker-Planck equations are identical (see the discussion of Section 13^ . 

Similarly, in the spin case we also choose a bilineal coupling, F = ^jSj, with 
some real constants. For spin systems [Hs, F] is in general ifg-dependent and a generic 
form of the time evolution equation cannot be given. For the sake of simplicity we take 
Hs = —BiSi, i.e. the isotropic spin {Dij = in (H2l)). This is the closest analogue to 
the open quantum oscillator discussed above (section 13^ . Then, the quantum master 
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equation (1461) follows in phase-space from 
1 d 



S dm 



m X B 



cff J 



m X A 



m X ( Boff - T 



d_ 
dm 



M X Beff 



W^j") (50) 



where [Befrjj = SBj is the classical effective field [see Eq. (fT9!) ]. A has components 
T 1 dF dF 

= -sd^,d^, ' ^''^ 

(if F is linear, then Kjk = S'j^j^j) and M = m{fj["'^ — S) + i(m x Ij)fj^\ Observe that 
M = 0(1/5*), so the last term vanishes as S* ^ oo recovering the classical Fokker-Planck 
equation [26] . 

Hence, even in the case where both system (Hg) and coupling (F) are linear in 
the spin variables, classical and quantum Fokker-Planck equations are different, unlike 
the open harmonic oscillator. Within the phase-space formalism this difference can be 
worked out explicitly. 



4-2. The classical limit of [J^l- general case 

To finish, let us go back to (l49ll and compute its classical limit for general Hs and F. 
Using ( !39l) one finds that. 



Im' 



5^ 1 ^ dF dF ^ 
AS"^ dnij drrik ^ 



1^{\H,[S^-^),F{S^^^)]) 



-1/S{H,,F} 



Hence, 



Im(F(5('^)))lm([/Js(5('^)),F(5('^))]) ^^^l^i^-F} 



Since {Hs,F} = ^eijk{dm,Hs){dmkF)mi we finally obtain: 



1 d 
S dm 



[m X Beff) — m X A 



m X ( B 



>cfr 



T 



d 
dm 



e 



(52) 
(53) 

(54) 
(55) 



where Ajk has been defined in ( 1511) and Beg in (]T9l) . Eq. (1551) is nothing but the classical 
Fokker-Planck equation for classical spins [cf. Eq. fl50|) ] \37\ (see also [26]). 



5. Conclusions 



The phase-space formulation of Quantum Mechanics provides deep insight into the 
quantum-classical correspondence, mainly due to the fact that the mathematical nature 
of the observables does not change when going from the classical to the quantum theory: 
they are always functions on the phase space. 

It is well-known that for non-relativistic spinless particles Bopp operators simplify 
the manipulations in the phase space formalism [20l [7] . In this paper we have generalized 
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Bopp operators for spin systems and applied them to the problem of the classical limit 
of spins as well as to open quantum spin systems. 

The dynamics of open quantum systems has received a renewed attention due to 
its role in explaining the emergence of the classical world from quantum mechanics |38j . 
Open quantum systems are handled with equations for the reduced density matrix, cf. 
Eq. fH6l) . They consist of the unitary evolution (von Neumann) and a non-unitary 
term which incorporates the bath influence. For the spinless case the transformation 
of the time evolution equations into phase space yields the quantum generalization of 
the corresponding Fokker-Planck equations. We have accomplished this task in the spin 
case with the help of the generalized Bopp operators. We have worked out the classical 
limit {S —>■ oo) of the quantum master equation for the density matrix, and recovered 
the classical Fokker-Planck equations. 
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Appendix A. 



We present here the details for the calculation of Eq. (1251) . For that we compute 

S3 



exphcitely Wg''^ ^Yim, using fl2Tl) . We first notice that [T9] : 



W, 



S3 



(A.l) 



Using that F{C?) cosO = F{1) cos9 and that in the sum of fl2Tl) j^ax = 1, together with 

1 ao 



ao 



(25+1)! 



ai 



we can write. 



(0 



COS 6Y1. 



Im 



25 + 2 



25 + 2 



5+« cos0U5-«n 



'■Im 



here S^^^^ acts only over cos 6^ and S over Yim with, see Eq. (I24 



86 sin 6 
Next we make use of the relations: 



cos OYim = aiYl+i^rn + a2Yi_i^, 

d 
o9 



(A.2) 



(A.3) 



(A.4) 



(A.5) 



Bopp operators and phase-space spin dynamics 
with 



(/ - m + + m + 1) 
(2/+ l)(2/ + 3) 



^2 



(/ — m)(/ + m) 
(2/ - 1)(2/ + 1) 



/3i = ai/, (32 = -02(1 + I). 
Manipulating the terms entering (1A.3P we find: 



F^-''{l)F''-^ 



cos 9Y, 



Im 



= a,F'-''{l)F'^-\l + l)Yi+,^ 
+ a2F'-''{l)F"-\l-l)Yi^r, 



Im 



and 



Now, taking into account that 

1 



we get: 
and 



Yl-lm 

COS OYi. 



- (32F^~''{l)F''-\l-l)Yi^ 



d 



p^^((i+l)cos(- + sin«|)y,„ 



Im 



= F^-"(/)F"-^(/ + l) 
+ F^-"{l)F''-\l - 1) 



2/ + 1 

1 

2/ + 1 



9 

(/ + l)cos^ + sin^— 

O'c' 



1^ 



Im 



d 

IcosO — sin^— 
06 



Y, 



Im 



cos^U5-«r^™ 



- F^-''(l)F''-\l + l) 
+ 1) 

+ 177- >«m- 



2/ + 1 
/+ 1 

21 + 1 



9 

(/ + l)cos^ + sin0— 

Ov 

d 

(/ + l)cos^ + sin0— 
o't7 



Ylm 



Finally, we obtain: 



with r^f \ i = l,2 defined in dZ?]). 



06 



\_d_ 

2d^. 



Y 



Im 
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